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Objectives 


♦ Determine the suitability of various inviscid flux schemes for 
DNS/LES of supersonic turbulent boundary layers - foundational for 
more complex turbulent flows in aerospace applications: 

• Shock-turbulent boundary layer interactions 

• Supersonic film cooling 




2-D Expansion-Compression Corner 



2-D Supersonic Film Cooling 
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Approach 


♦ A number of approaches for DNS or LES of supersonic turbulent 
boundary layers have been used over the past two decades 

• Fully spatially evolving simulations 

- Forced transition from laminar inflow BC 

- Quasi-turbulent inflow BC 

- Recycling of inflow BC from point further downstream 

• Extended temporally-developing simulations 

• Temporally-developing simulations 

♦ Temporally-developing approach used here 

• Inexpensive: relatively small domain can be used 

• Allows spatial averaging and periodic boundary conditions in both 
streamwise (x) and spanwise (z) directions 

• Significant disadvantage is that the flow field is never truly statistically 
stationary - rather, the boundary layer slowly grows with time. 
However, Martin (2004, 2007) and Xu and Martin (2004) have shown 
that if a relatively short time window is chosen, good results can be 
obtained. 


Temporally-Developing Turbulent Boundary Layer 
at Mach 2.889, Re = 1 .93 x 1 0 6 nr 1 


Turbulent Boundary Layer Flow - Synthetic Schlieren Image 



♦ t* is a normalized timescale t IL / 9.6L (domain flowthrough timescale) 


Scope 


♦ Emphasis in this paper is on studying the basic resolving 
characteristics of spatial discretizations and inviscid flux 
schemes for this compressible flow problem, not on LES 
sub-grid scale (SGS) models 

♦ Grid resolutions considered here are from DNS levels to fine- 
grid LES levels 

♦ Adiabatic wall BC selected as this case is well-understood 
and Morkovin’s hypothesis holds at Mach 5 and below 

♦ Viscous discretization is 2 nd -order accurate only (limited 
tests with higher-order accuracy do not appear to change 
results appreciably) 

♦ Time advancement is 4-stage Runge-Kutta, limited by 
stability restrictions to time steps significantly smaller than 
those needed to resolve turbulence 
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Governing Equations of Fluid Mechanics 


Navier-Stokes Equations in Terms of Conservative Variables: 
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Conservative Differencing Forms the Flux at a Half-Point: 


dU dF(U) 

— + = 0 

dt dx 



dF(U) K 

— I=|> AU, = -(At/Ax)(F i+y , - F„ /2 ) 


6 



Skew-Symmetric Central Difference Schemes 
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SSCD-2 F j+ y 2 - F avg (i,i+1 ) 

SSCD-4 F i+1/2 = (4/3)F avg (i,i+1) - (1/6)[F avg (i-1 j+1) + F avg (i,i+2)] 

SSCD-6 F i+1/2 = (3/2)F avg (i,i+1 ) - (3/1 0)[F avg (i-1 j+1 ) + F avg (i,i+2)] 

+(1/30)[F avg (i-2,i+1) + F avg (i-1 ,i+2) + F avg (i,i+3)] 

SSCD-8 F i+1/2 = (1 6/1 0)F avg (i,i+1 ) - (4/1 0)[F avg (i-1 j+1 ) + F avg (i,i+2)] 

+(8/1 05)[F avg (i-2,i+1 ) + F avg (i-1 ,i+2) + F avg (i,i+3)] 
-(1/140)[F avg (i-3,i+1) + F avg (i-2,i+2) + F avg (i-1 ,i+3) + F avg (i,i+4)] 


^avg(^ — Pavg ^avg ^avg ^avg 

where V avg = [1 , U av g, V av g, w avg , © avg + /2(u avg + V avg + w avg ) + Pavg^Pavg] 

and P avg = [0, P avg , 0, 0, 0] 


For any variable cp, cp avg = 1 / 2 (cp i1 + cp i2 ). This F avg is very similar to the skew- 
symmetric scheme of Kennedy and Gruber (2008) (Also see Pirozzoli (2010)). Note 
if F avg = (F n + F i2 ), then standard divergence form central differencing formulas result. 



Filtered Central Difference Scheme 


Eighth-Order Divergence Form Central Difference Scheme 

DCD-8 F i+1/2 = [-3Fj_ 3 + 29Fj. 2 - 1 39F M + 533Fj + 533F i+1 - 1 39F i+2 +29F i+3 - 3F i+4 ]/840 


Combined with Compact Filtering Algorithm 
(Visbal and Gaitonde, 1999) at interior points, 
and High Order Boundary Filtering Algorithm 
(Gaitonde and Visbal, 2000) at the lower and 
upper domain boundaries. 

The filter completely filters wavelengths at the 
spatial Nyquist frequency k = tt (2 p.p.w.) 

The filter transfer function is > 0.999 at k = tt/2 
(4 p.p.w.), and is > 0.99 at k = 2 tt/ 3 (3 p.p.w.) 

Filter is applied after one full Runge-Kutta 
timestep 


Filter Transfer Function at Interior Points 



Combination of DCD-8 with filtering algorithm 
is referred to DCD-8-F in this work 
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Upwind-Biased Roe Schemes 
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Roe Flux using “left” and “right” interpolations for the half-node 

Fj+i/2 = y^L + F i+1/2 , R ) - V 2 1 A roe | ( U j + y 2j R - Uj +1/2jL ) 
where F i+1/2jL and U i+y2jL are formed from [p, u, v, w, P] i+ i /2jL 
and F j+1/2R and U i+ y 2jR are formed from [p, u, v, w, P] i+1/2 R 

Primitive variable reconstruction for generic flow variable q> (no slope limiting) 

UBR-3 cp i+1/2L = (-cp M + 5cp, + 2cp i+1 )/6 

UBR-5 cp i+1/2L = (4cpj. 2 - 26<p M + 94cpj + 54cp i+1 - 6cp i+2 )/1 20 

UBR-7 cp i+ i/ 2jL = (-6cpj_ 3 + 50cpi_ 2 - 202cpj_-, + 638cpj + 428cp i+1 - 76cp j+2 + 8cp j+3 )/840 

The cp i+ y 2 , R values are formed from a flipped interpolation 

Note: none of the inviscid flux schemes used here, including the upwind- 
biased methods, should be considered as “shock-capturing” (i.e. MUSCL or 
WENO methods) 


Modified Scaled Wavenumber, 


Fourier Analysis Characteristics 


Real Part of Modified Scaled Wave- Imaginary Part of Modified Scaled 

number: Dispersion (Phase) Wave-number: Dissipation 




♦ The dispersion (phase) characteristics of each UBR scheme are the same as the SSCD/DCD 
scheme of one higher order of accuracy 

♦ SSCD-2 has a limited range of phase accuracy, higher-order schemes get progressively better 

♦ Central schemes are non-dissipative, while all upwind schemes experience increasing 
dissipation error at higher wavenumbers 

♦ Higher-order upwind schemes have smaller dissipation error at lower wavenumbers than lower- 
order upwind schemes 
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Additional Numerical Details 


Air at T„ = 298.15 K, R = 287 J/Kg-K, y = 1.4, c p = 1004.5 J/Kg-K, 
a TO = 346.117 m/s, p = Sutherland relation, A = Sutherland relation 

Domain: 0 < x < 9.6L, 0 < y < 10L, 0 < z < 4.8L, L = 0.01 m 

Grid A: 128x97x128 
Grid B: 192x129x192 
Grid C: 256x145x256 

uniform spacing and periodic B.Cs. in x- and z-directions 

nonuniform spacing and adiabatic viscous wall B.C. (lower boundary) and 
non-reflecting characteristic B.C. (upper boundary) in y-direction 

4-stage Runge-Kutta time advancement 

Initial condition obtained from prescribed turbulent mean flow profile and 
velocity fluctuations obtained from channel flow simulations (see paper) 


ummary of Freestream Conditions and Cases Run 
Nominal Reynolds Number Cases 



M w 

P«, 

Re« 


SSCD 



UBR 


DCD-8-F 

(m/s) 


(Pa) 

(nr 1 ) 

-2 

-4 

-6 

-8 

-3 

-5 

-7 


50 

0.144 

25331 

8.05 xiO 5 

B 

B 

B 

B 

B 

B 

B 

B 

250 

0.722 

5066 

8.05 x 10 s 

B 

B 

B 

B 

B 

B 

B 

B 








A 



A 


500 

1.445 

3800 

1.21 x 10 6 

B 

B 

B 

B 

B 

B 

B 

B 








C 



C 









A 



A 


750 

2.167 

3800 

1.81 xiO 6 

B 

B 

B 

B 

B 

B 

B 

B 








C 



C 









A 



A 


1000 

2.889 

4560 

2.90 xiO 6 

B 

B 

B 

B 

B 

B 

B 

B 


c c 


Grid A = 1 28 x 97 x 1 28 Ay w + = 0.33, Ax + = 24, Az + = 1 2 
Grid B = 1 92 x 1 29 x 1 92 Ay w + « 0.33, Ax + * 1 6, Az + * 8 
Grid C = 256 x 1 45 x 256 Ay w + « 0.33, Ax + * 1 2, Az + * 6 


Note: only 
SSCD-8 and 
UBR-7 were run 
on grids A and C 
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ummary of Freestream Conditions and Cases Run 
Nominal Reynolds Numbers Cases (Cont’d) 



M„ 

Poo 

Re« 

SSCD 



UBR 


DCD-8-F 

(m/s) 


(Pa) 

(nr 1 ) 

-2 -4 -6 

-8 

-3 

-5 

-7 










A 


1250 

3.612 

5066 

4.03 xio 6 

B 

B 

B 

B 

B 

B 






C 



C 










A 


1500 

4.334 

6333 

6.04 xio 6 



B 

B 

B 

B 









C 


1750 

5.056 

7599 

8.46 xio 6 



B 

B 

B 

B 


Note: only SSCD-8 and UBR-7 were run on grids A and C 

Grid A = 1 28 x 97 x 1 28 Ay w + * 0.33, Ax + * 24, Az + * 1 2 

Grid B = 1 92 x 1 29 x 1 92 Ay w + * 0.33, Ax + * 1 6, Az + * 8 

Grid C = 256 x 1 45 x 256 Ay w + * 0.33, Ax + * 1 2, Az + * 6 



ummary of Freestream Conditions and Cases Run 
Reduced Reynolds Number Cases 



M„ 

Poo 

Re„ 

SSCD 


UBR DCD-8-F 

(m/s) 


(Pa) 

(nr 1 ) 

-2 -4 -6 

-8 

-3 -5 -7 






A 

A 

500 

1.445 

2533 

8.05 xIO 5 


B 

B 






C 

C 






A 

A 

750 

2.167 

2533 

1.21 xiO 6 


B 

B 






C 

C 






A 

A 

1000 

2.889 

3040 

1 .93 x 1 0 6 


B 

B 


Note: only SSCD-8 and UBR-7 were run at these conditions 

Grid A = 1 28 x 97 x 1 28 Ay w + * 0.23, Ax + * 1 6, Az + * 8 
Grid B = 1 92 x 1 29 x 1 92 Ay w + * 0.22, Ax + *11, Az + « 5 
Grid C = 256 x 1 45 x 256 Ay w + * 0.22, Ax + * 8, Az + * 4 



ummary of Freestream Conditions and Cases Run 
Reduced Reynolds Number Cases (Cont’d) 



M„ 

Poo 

Re« 

SSCD 


UBR DCD-8-F 

(m/s) 


(Pa) 

(nr 1 ) 

-2 -4 -6 

-8 

-3 -5 -7 






A 

A 

1250 

3.612 

3378 

2.68 xio 6 


B 

B 






C 

C 







A 

1500 

4.334 

4222 

4.03 x 10 6 


B 

B 






C 

C 







A 

1750 

5.056 

5066 

5.64 x 1 0 6 



B 


Note: only SSCD-8 and UBR-7 were run at these conditions 


Grid A = 1 28 x 97 x 1 28 Ay w + * 0.23, Ax + * 1 6, Az + * 8 
Grid B = 1 92 x 1 29 x 1 92 Ay w + * 0.22, Ax + *11, Az + « 5 
Grid C = 256 x 1 45 x 256 Ay w + * 0.22, Ax + * 8, Az + * 4 
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urbulent Boundary Layer at Mach 2.889, Re = 1.93 x 10 6 nr 1 

Iso-surfaces of Q-criterion 


Normalized Velocity 


Normalized Temperature 




Turbulent Boundary 


Layer Flow - Iso-surfaces of Q = 0.1 (U x /8 99 )' 


Y 


Normalized Time: 1 1 .50 


2.50 

2.10 

1.70 

1.30 

0.90 


Normalized Time: 1 1 .50 


♦ Simulation run using SSCD-8 scheme on Grid C: 256 x 145 x 256 points 

♦ Iso-surfaces colored by normalized velocity magnitude U/IL and temperature T/T 

♦ t* is a normalized timescale t U M / 9.6L (domain flowthrough timescale) 

♦ Turbulent flow structure is evident, as are sharp gradients near wall 
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urbulent Boundary Layer at Mach 2.889, Re = 1.93 x 10 6 nr 
Contour Plots in x-z plane at y + = 12.6 


i 


Normalized Velocity 


Normalized Temperature 


Turbulent Boundary 


Layer Flow - Velocity Co itours at y + = 12.6 


Turbulent Boundary 


Layer Flow - Temperature Contours at y + = 12.6 


z 



X 




Normalized Time: 1 1 .50 



♦ Velocity magnitude and temperature are normalized by freestream values 

♦ This x-z cutting lane at y + =12.6 is the location of peak velocity fluctuations 

♦ Elongated streaky structures characteristic of near-wall turbulence are easily 
observed, along with sharp fluctuations in temperature 
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Normalized Streamwise Velocity, u/U 



urbulent Boundary Layer at Mach 2.889, Re = 1.93 x 10 6 nr 

Scatter Plots in x-z plane at y + = 12.6 
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Normalized Streamwise Velocity 



z/L 


Normalized Temperature 
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♦ Streamwise velocity and temperature are normalized by freestream values 

♦ The green line and symbols indicates a results extracted from a single line at x / L = 4.8 

♦ The blue lines indicate the average value as well as ± one standard deviation 

♦ The streaky structures are quite evident in the constant x-station spanwise profiles 

♦ As the calculation of the inviscid flux involves computing derivatives of variables 
involving p, u, v, w and p, the need for adequate grid resolution and high numerical 
accuracy in computing these types of flows is obvious 


Variable Normalized by Average Value 


urbulent Boundary Layer at Mach 2.889, Re = 1.93 x 10 6 nr 1 
Normalized Fluctuations along X / L = 4.8 Line at y + = 12.6 


Fluctuating u and T 



Fluctuating Mach Numbers 



♦ In left panel, normalized streamwise velocity and temperature fluctuations are 
compared spanwise along the x / L = 4.8 line at y + = 12.6 

♦ The two fluctuations are generally anti-correlated. Also shown is a plot of the 
“fluctuation” form of the Strong Reynolds Analogy for comparison. At the y + = 12.6 
plane, the value of the RMS form of the SRA is 0.897 - a strong correlation on average. 

♦ In the right panel, two expressions for the fluctuating Mach number are compared. 

The effect of the temperature/sound speed fluctuations (anti-correlated with u) is 
evident. 


Skin Friction Coefficient, c f 


Comparison of Non-Dissipative SSCD Schemes 
at Mach 0.144, Re = 8.05 x 10 5 nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



II) Karman-Schoenherr relation for skin friction at this condition 

♦ All four similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile of Simens et al. (2009) and Jimenez et al. (2010) at Re 52 = 1551 
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Skin Friction Coefficient, c f 


Comparison of Non-Dissipative SSCD Schemes 
at Mach 0.722, Re = 8.05 x 10 5 nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



♦ All four SSCD schemes generally compare quite well with the transformed (van Driest 
II) Karman-Schoenherr relation for skin friction at this condition 


♦ All four similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551 
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Skin Friction Coefficient, c f 


Comparison of Non-Dissipative SSCD Schemes 
at Mach 1.445, Re = 1.21 x 10 6 nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



Driest II) Karman-Schoenherr relation for skin friction at this condition. SSCD-2 
exhibits a skin friction slightly higher than the relation. 

♦ All three high-order methods similarly match the laminar sublayer, log law, and 
incompressible DNS velocity profile at Re 52 = 1551. SSCD-2 begins to exhibit a slightly 
smaller wake strength than the others. 
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Skin Friction Coefficient, c f 


Comparison of Non-Dissipative SSCD Schemes 
at Mach 2.167, Re = 1.81 x 10 6 nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



♦ The SSCD-4, -6 and -8 schemes generally compare quite well with the transformed (van 
Driest II) Karman-Schoenherr relation for skin friction at this condition. SSCD-2 
exhibits a skin friction higher than the relation. 


♦ All three high-order methods similarly match the laminar sublayer, log law, and 
incompressible DNS velocity profile at Re 52 = 1551. SSCD-2 exhibits a smaller wake 
strength than the others, and also starts to fall slightly below the log law line. 
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Skin Friction Coefficient, c f 


Comparison of Non-Dissipative SSCD Schemes 
at Mach 2.889, Re = 2.90 x 10 6 nr 1 on Grid B 


Skin Friction Mean Velocity Profiles 



♦ The SSCD-4, -6 and -8 schemes generally compare quite well with the transformed (van 
Driest II) Karman-Schoenherr relation for skin friction at this condition. SSCD-2 
exhibits a skin friction higher than the relation. 


♦ All three high-order methods similarly match the laminar sublayer, log law, and 
incompressible DNS velocity profile at Re 52 = 1551. SSCD-2 exhibits a smaller wake 
strength than the others, and also starts to fall slightly below the log law line. 
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Skin Friction Coefficient, c f 


Comparison of Non-Dissipative SSCD Schemes 
at Mach 3.612, Re = 4.03 x 10 6 nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



♦ The SSCD-6 and -8 schemes generally compare quite well with the transformed (van 
Driest II) Karman-Schoenherr relation for skin friction at this condition. SSCD-2 and -4 
were not stable at this condition. 


♦ The two high-order methods similarly match the laminar sublayer, log law, and 
incompressible DNS velocity profile at Re 52 = 1551 

♦ Note that the boundary layer thickness (at equivalent Re 52 ) has grown noticably 
shorter in wall coordinates than was the case for the subsonic cases 



Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 0.144, Re = 8.05 x 10 s nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



(van Driest II) Karman-Schoenherr relation for skin friction at this condition 

♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551 

♦ UBR-5 does less well, and UBR-3 does poorly. This can be traced to the excessive 
dissipation in the basic Roe scheme at low Mach numbers. Fixes for this published in 
recent years bring UBR-3 close to UBR-5 for this case. 


Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 0.722, Re = 8.05 x 10 s nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



♦ The UBR-7 and DCD-8-F schemes generally compare quite well with the transformed 
(van Driest II) Karman-Schoenherr relation for skin friction at this condition 

♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551 

♦ UBR-5 is now nearly as good as UBR-7 and DCD-8-F, rising slightly over the log law line 

♦ UBR-3 is better than it was at Mach 0.144, but still has an unacceptable velocity profile 
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Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 1.445, Re = 1.21 x 10 6 nr 1 on Grid B 


Skin Friction Mean Velocity Profiles 



(van Driest II) Karman-Schoenherr relation for skin friction at this condition 

♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551. Both do not predict quite the same wake strength that SSCD-8 
did. 

♦ UBR-5 is nearly as good as UBR-7 and DCD-8-F 

♦ UBR-3 is better than it was at Mach 0.144, but still has an unacceptable velocity profile 


Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 2.167, Re = 1.81 x 10 6 nr 1 on Grid B 


Skin Friction Mean Velocity Profiles 



♦ The UBR-7 and DCD-8-F schemes generally compare well with the transformed (van 
Driest II) Karman-Schoenherr relation for skin friction at this condition 


♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551. Both do not predict quite the same wake strength that SSCD-8 
did. 


♦ UBR-5 is nearly as good as UBR-7 and DCD-8-F 

♦ UBR-3 is better than it was at Mach 0.144, but still has an unacceptable velocity profile 


Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 2.889, Re = 2.90 x 10 6 nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



♦ The UBR-7 and DCD-8-F schemes generally compare quite well with the transformed 
(van Driest II) Karman-Schoenherr relation for skin friction at this condition 


♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551. UBR-7 does not predict quite the same wake strength that 
SSCD-8 did. 

♦ UBR-5 is nearly as good as UBR-7 and DCD-8-F 

♦ UBR-3 is better than it was at Mach 0.144, but still has an unacceptable velocity profile 


Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 3.612, Re = 4.03 x 10 6 nr 1 on Grid B 

Skin Friction Mean Velocity Profiles 



♦ The UBR-7 and DCD-8-F schemes generally compare quite well with the transformed 
(van Driest II) Karman-Schoenherr relation for skin friction at this condition 

♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551. UBR-7 does not predict quite the same wake strength that 
SSCD-8 did. 

♦ UBR-5 is nearly as good as UBR-7 and DCD-8-F 

♦ UBR-3 is better than it was at Mach 0.144, but still has an unacceptable velocity profile 


Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 4.334, Re = 6.04 x 10 6 nr 1 on Grid B 


Skin Friction Mean Velocity Profiles 



♦ The UBR-7 and DCD-8-F schemes generally compare quite well with the transformed 
(van Driest II) Karman-Schoenherr relation for skin friction at this condition 


♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551. DCD-8-F exhibits a slightly higher wake strength than UBR-7. 

♦ UBR-5 is nearly as good as UBR-7 and DCD-8-F 

♦ UBR-3 is better than it was at Mach 0.144, but still has an unacceptable velocity profile 
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Skin Friction Coefficient, c f 



Comparison of Dissipative UBR and DCD-8-F Schemes 
at Mach 5.056, Re = 8.46 x 10 6 nr 1 on Grid B 


Skin Friction Mean Velocity Profiles 



♦ The UBR-7 and DCD-8-F schemes generally compare quite well with the transformed 
(van Driest II) Karman-Schoenherr relation for skin friction at this condition 


♦ They similarly match the laminar sublayer, log law, and incompressible DNS velocity 
profile at Re 52 = 1551. DCD-8-F exhibits a slightly higher wake strength than UBR-7. 

♦ UBR-5 is nearly as good as UBR-7 and DCD-8-F 


♦ UBR-3 is better than it was at Mach 0.144, but still has an unacceptable velocity profile 
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Fluctuating RMS Mach Number, VM’ 


Fluctuating RMS Mach Number Profiles 
Nominal Reynolds Number Cases on Grid B 


SSCD-8 



Distance from Wall, y/5 99 



Distance from Wall, y/5 99 


♦ Clear increase in fluctuating Mach number with freestream Mach number 

♦ Some differences in the profiles are evident between the two numerical methods at 
supersonic conditions. UBR-7 predicts a higher fluctuating Mach number near the wall 
than SSCD-8 does. 

♦ The results suggest that fluctuating Mach numbers above 0.34 are problematic for 
non-dissipative schemes, while dissipative schemes can handle values up to 0.5 
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Skin Friction Coefficient, c f 


Effect of Grid Resolution on Skin Friction and 
Velocity Profiles at Mach 3.612 


Skin Friction Mean Velocity Profiles 



♦ Grid C at reduced Reynolds number had Ay w + ~ 0.22, Ax + ~ 8, Az + ~ 4 

♦ Little difference between the two grid resolutions (in wall coordinates) is evident in the 
results for SSCD-8 

♦ Higher grid resolution (in wall coordinates) does appear to improve the prediction of 
skin friction for UBR-7, and also increase the wake strength 
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Van Driest Scaled RMS Velociffes 


Effect of Grid Resolution on 
RMS Velocity Components at Mach 3.612 


Comparison in Wall Coordinates 


Comparison in Outer Coordinates 



1.0 - 
0.5 - 


10 ° 10 1 _ + 10 2 
Distance from Wall, y 


4.0 r 
3.5 - 

3.0 - 


2.5 - 
2.0 - 


Incompressible DNS, Re 0 = 1551 
SSCD-8, Grid B, Re = 4.03 x 10 6 , Re 62 = 1553 
SSCD-8, Grid C, Re = 2.68 x 10 6 , Re* = 1553 
UBR-7, Grid B, Re = 4.03 x 10 6 , Re„ = 1549 


UBR-7, Grid C, Re 


. ..+ 



♦ Grid B at nominal Reynolds number had Ay w + ~ 0.33, Ax + ~ 16, Az + » 8 

♦ Grid C at reduced Reynolds number had Ay w + ~ 0.22, Ax + » 8, Az + » 4 

♦ UBR-7 at either grid resolution (in wall coordinates) appears to predict higher near-wall 
peak values of u vd ,+ than SSCD-8 does - the difference may narrow at higher resolution 

♦ Note comparatively lower values of u vd ,+ in some of the outer layer (wake region) at the 
coarser grid resolution - may possibly help explain reduced wake strength 
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Coles' Wake Strength Parameter, 



Effect of Grid Resolution on 
Coles’ Wake Strength Parameter 


Nominal Reynolds Number Cases 
(Ay w + * 0.33, Ax + * 16, Az + * 8) 



Reduced Reynolds Number Cases 
(Ay w + * 0.22, Ax + * 8, Az + * 4) 



♦ The wake strength is the difference between the velocity profile and the log law line at 
the outer edge of the boundary layer - here the 99% boundary layer thickness was used 


♦ Coles’ Profile for wake strength as a function of Re 0 was compiled from a large set of 
the available experimental data in the early 1960s 

♦ The differences between the SSCD-8 and UBR-7 schemes at the coarser grid resolution 
are evident. The differences narrow considerably at higher grid resolution. 
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Shape Factor, 


Effect of Grid Resolution on 
Boundary Layer Shape Factor 


Nominal Reynolds Number Cases Reduced Reynolds Number Cases 

(Ay w + * 0.33, Ax + * 16, Az + * 8) (Ay w + * 0.22, Ax + * 8, Az + * 4) 



♦ The boundary layer shape factor, H, is the ratio between the displacement thickness, 5*, 
and the momentum thickness, O 


♦ Michel’s curve for H (Gatski, 2009) is simple for an adiabatic turbulent boundary layer: 
H = H inc + 0.4 MJ 

♦ Both numerical schemes show excellent agreement with Michel’s curve. Very little 
effect from grid resolution is evident. 
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Effect of Assuming Constant c p with 
Sutherland Formulas for p and A 



♦ This combination of assumptions leads to unrealistically low Prandtl numbers at 
elevated temperatures - should have used variable c p , or assumed constant Pr and 
calculated A from it (A = p c p / Pr) 
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Recovery Factor, r 


Effect of Grid Resolution, Numerical Scheme and A on 
Turbulent Boundary Layer Recovery Factor 


SSCD-8 


UBR-7 
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♦ The turbulent boundary layer recovery factor is a measure of how much of the total 
temperature of the flow is actually recovered at the wall in adiabatic flow 

♦ An additional 11 simulations were run at the highest grid resolution with the 
assumption of fixed Pr = Pr (298.15 K), and A = p c p / Pr 

♦ The effect of this assumption on recovery factor is significant 

♦ Finer wall-normal grid spacing appears to help bring the observed recovery 
temperature closer to expected values 

♦ Also note that UBR-7 is typically slightly better here than SSCD-8 
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Summary 


♦ Four non-dissipative skew-symmetric central difference (SSCD-2, SSCD-4, SSCD-6 
and SSCD-8) schemes, three upwind-biased Roe (UBR-3, UBR-5 and UBR-7) 
schemes, and one filtered divergence form central difference (DCD-8-F) scheme 
were evaluated on a temporally-developing compressible turbulent boundary layer 
problem. 

♦ The schemes were tested at Mach numbers ranging from effectively incompressible 
(Mach 0.144) to hypersonic (Mach 5.056). Three different grid sizes were used, and 
case sets at two different unit Reynolds numbers were run for all of the supersonic 
cases. The schemes were evaluated in terms of agreement with well-established 
skin friction and velocity profiles. 

♦ There are clear benefits to higher-order methods for this problem. Among the non- 
dissipative methods, at the baseline grid resolution in wall coordinates used for 
much of this work (Ax + » 16, Az + ~ 8), the SSCD-2 scheme clearly became less 
accurate as the Mach number increased. Further, both SSCD-2 and SSCD-4 became 
unstable for cases at Mach 3.612 and higher. SSCD-6 and SSCD-8 remained stable 
at Mach 3.612, and SSCD-8 was stable at Mach 4.334 at the two highest grid 
resolutions (in wall coordinates) tested. The results suggest that fluctuating Mach 
numbers above 0.34 become problematic for a non-dissipative scheme, though it 
appears that fluctuating Mach numbers of 0.43 can be tolerated with higher grid 
resolution. 
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Summary 


♦ Among the dissipative schemes, UBR-7 offers clear improvements over UBR-3, and 
to a lesser extent UBR-5, across the entire Mach number range. UBR-3 performed 
particularly poorly at the lowest Mach number tested, likely due to the relatively 
higher dissipation characteristics of the Roe solver at low speeds. The DCD-8-F 
scheme appeared to perform similarly to UBR-7, and possibly slightly better. UBR-7 
and DCD-8-F appear to have an excellent combination of accuracy and robustness. 
The dissipative UBR-7 scheme handled fluctuating Mach numbers up to 0.5 without 
difficulty. 

♦ The SSCD-8 and UBR-7 schemes generally performed very similarly on this problem. 
Some differences between these schemes were noted with respect to turbulence 
intensities and wake strength at nominal grid resolution in wall coordinates (Ax + ~ 8, 
Az + ~ 4). These differences diminished considerably at higher grid resolution (Ax + ~ 
8, Az + » 4). The SSCD-8 scheme at both grid resolutions, and the UBR-7 scheme at 
higher grid resolution, exhibited reasonably good agreement with Coles' fit for wake 
strength as a function of Re 0 . The SSCD-8 and UBR-7 results also agree well with 
Michel's formula for boundary layer shape factor as a function of Mach number. 
Some small differences between the two schemes were noted for the turbulent 
boundary layer recovery factor. A reduction in Ay w + from » 0.33 to » 0.22 also had a 
significant effect on the recovery factor. 
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DNS 


LES/ILES 


HRANS-LES 


RANS 


Turbulence Modeling Strategies 


Directly resolve all scales of turbulent motion. There is no 
eddy viscosity, and the intrinsic dissipation of the inviscid 
flux must be « laminar viscosity. 


Directly resolve largest scales of turbulent motion, and 
model the smaller scales using one of two approaches: 

A) classic LES - a subgrid-scale (SGS) turbulence model, 
requiring that the intrinsic dissipation of the inviscid flux is 
« SGS eddy viscosity, or B) ILES - where the intrinsic 
dissipation of a 2 nd order accurate inviscid flux 
approximately mimics the SGS eddy viscosity. 


Attempts to directly resolve turbulence only in regions 
with adequate grid resolution, otherwise turbulence is 
modeled. In LES regions the algorithmic requirements 
would be consistent with the above description. 


Turbulence is entirely modeled. 
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